Spatial solitons and modulational instability in the presence of large birefringence: 
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With reference to spatially non-local nematic liquid crystals, we develop a theory of optical spa- 
tial solitons and modulational instability in anisotropic media with arbitrarily large birefringence. 
Asymmetric spatial profiles and multivalued features are predicted for self-localized light versus 
walk-off angle. The results hold valid for generic self-focusing birefringent media and apply to large 
angle steering of individual and multiple self-trapped optical beams. 



An appealing approach towards the realization of dig- 
ital multidimensional all-optical processors and informa- 
tion routers is the use of spatially self-trapped optical 
filaments -or solitons- as readdressable light pencils able 
to guide signals in arbitrary directions. Although recent 
years have witnessed widespread investigations of spa- 
tial solitons in various nonlinear systems Q, j t date a 
proper description and, hence, modeling and prediction 
of such functionality is hampered by two main issues. 
The first is somewhat of a technical origin: while the most 
studied nonlincarities for stable two-dimensional spatial 
solitons rely on birefringent materials, Q the effects of 
anisotropy have been accounted only for small walk-off or 
propagation along one of the principal axes (e.g. in pho- 
torefractives or quadratic media). Dealing 
with anisotropy in a perturbative way, current models fail 
to predict the formation of self-collimated beams read- 
dressable over wide-angles. The second issue relates to 
the ubiquitous paraxial approximation, which is unable 
to treat propagation at large angles with respect to the 
input wavevector, i.e. to the launch direction z of the 
beam generating the soliton. 

The above considerations hold valid also for spatial 
modulational instability (MI) , a process tipically accom- 
panying (or precurring) solitons. MI describes unstable 
plane-waves which, through self-focusing, break up into 
transversely periodic patterns eventually evolving into fil- 
aments. [lOj To date, the analysis of optical MI has been 
limited to small birefringence and paraxial behavior even 
in crystals with significant anisotropy. 

In this Letter we address nematic liquid crystals (NLC) 
as a natural environment to assess the role of a strong 
anisotropy in beam self-localization and MI. In doing so, 
by letting the beam be paraxial not with respect to z 
but to a rotated reference system, we develop a model 
which properly accounts for arbitrarily large walk-off and 
birefringence. While this allows to deal with beam steer- 
ing over large angles, 0] it also enlightens unexplored 
features of solitons and MI for arbitrary crystal orien- 



tations. The beam, polarized as an extraordinary wave, 
gives up radial symmetry while acquiring an asymmet- 
ric transverse profile. The latter, depending on angle 
of propagation, reveals a nontrivial distribution of both 
longitudinal and transverse components. Moreover, dif- 
fering anisotropic solitons or MI patterns can be expected 
for a given walk-off. 

The approach we introduce hereby is general and can 
be applied to any anisotropic nonlinear medium. For il- 
lustration sake and in order to pin-point a physically rel- 
evant system, we explicitly refer to a voltage-biased glass 
cell containing a thick layer of planarly-anchored nematic 
liquid crystals. Such configuration, encompassing a sig- 
nificant and externally-adjustable degree of birefringence 
(and walk-off), has been previously exploited for various 
experiments with spatial solitons (nematicons) and MI. 
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In the present context we adopt a refer- 
ence system and notation as in Fig. 1. The starting point 
is the vectorial wave equation: VxVxE = /c§e-E, with 
(constant) dielectric tensor s. Looking for a plane-wave 
solution E = A exp(i/co"2 ; ) propagating along z, the lin- 
ear homogeneous algebraic system 

C(n) ■ A — [n 2 (zz - I) + e] ■ A = 0, (1) 

where I is the identity matrix, z the unit vector in the 
z— direction, and zz the dyadic tensor whose elements are 
ZiZj with with i,j = {x, y, z}, yields the allowed values 
for n, i.e. ordinary and extraordinary refractive indices. 
Considering a light beam propagating in the midplane of 
a much thicker cell, for NLC with director n (e.g. mean 
orientation of molecular major axes) in the (y, z) plane, 
as in figure the relative permittivity tensor is given by 
Eij = £±5ij + Aenirij, and the anisotropy Ae constant in 
the illuminated region. The ordinary (o-) wave is polar- 
ized along x and n 2 = e±, while the extraordinary (e-) 
wave belongs to the (y, z) plane and the resulting index 
is n e (6 ) 2 = 2ej_(£_L + AE)/[2£x + AE + Aecos(20o)]- The 
unit- vector associated to the e-wave is denoted t(6*o), and 
its normal defines the walk-off direction s(8q). S(8q) is the 



2 




FIG. 1: Adopted coordinate system: The grey ellipse is a 
sketch of a LC molecule, n is its director, s and t are the 
walk-off direction and its normal, respectively. O = 9q in the 
absence of optical excitation, 5(6q) is the walk-off angle. The 
inset illustrates a feasible experimental arrangement, with X, 
Y and Z = z the axes in the laboratory frame. The applied 
voltage determines the elevation angle 770 of the molecules, 
which at zero bias lie in the Y, Z plane with azimuth £. The 
beam propagates along s in the plane (n,z). 



walk-off angle, with tan(<5) = [Aesm(0 o )cos(8 o )]/[£_i_ + 
Aecos(6'o) 2 ]. We omit hereafter the dependence on 6 . 

The general plane-wave solution with wavevector par- 
allel to z is a combination of e- and o-waves. In the fol- 
lowing, we only consider the e-wave polarized along t and 
propagating along s; the ordinary wave, being orthogonal 
to n, at the lowest order of approximation does not affect 
NLC molecular re-orientation through dipolc-ficld inter- 
action, because of the existence of a threshold known as 
the optical Freedericksz transition. (l4Hl5| Hence, for o- 
waves of intensity well below the Freedericksz threshold, 
the e-wave is the leading term in the NLC reorientational 
nonlinear response. 

The optical field perturbs the dielectric tensor as £ — > 
e + ( 2 8e, with e a smallness parameter to be taken equal 
to 1 at the end of the derivation. The expansion is written 
as E = [tE e + eF e + e 2 G e + ...] exp(i/c n e ,2), with E e , 
F e and G e depending on multiple slow scales 
t n = e n t and s n = e n s (n — 1,2,...) in the reference 
system (x,t,s). At the order 0(e), it is 

k 2 C(n e ) ■ F e = ik Q n e [z x Vi x (E e t) + Vi x z x (E e i)] - 



Of, 



ik a n e t[— 2 cos(<5) 
ik n e s[sm(S) §^ + cos(<5) 91 



fc n e x[sin(5)|^]- 

dti J' 



(2) 

The solvability condition implies the rhs of to 
be orthogonal to the null space of C(n e ), given by 
t: dE e /ds\ = 0. For the first-order vectorial cor- 
rection F, writing C(n e ) in the (x, t, s) system pro- 
vides: F* = (in e $m5/k \ x )dE e /dxi, F* = 0, F e s = 
(in e cos 5/ko\ s )dE e /dti; being X x ^ s the non-vanishing 
eigenvalues of C{n e ) (see Eqs. {J£J below). Hence, at this 
order of approximation, the electric field is not linearly 
polarized as an extraordinary wave, but its polarization 



varies across the finite trasverse profile. At the order 
0(e 2 ) 

k 2 £(n e ) ■ G e = -k 2 Q 8e ■ tE e + Vi x Vi x (E e i) 
ikon e [z x Vi x F e + Vi x z x F e ]+ (3) 
ik n e [z x V2 x (E e t) + V2 x z x (E e t)]. 

Using the result obtained at the previous order, from the 
solvability condition t • C(n e ) ■ G e = it is found (in the 
original scales, e — > 1) 



2lkoneCOsi§) — +Dt — r+Dx — 



k£(t-5e-t)E e = 
(4) 

i.e., the paraxial propagation equation in the walk-off 
system. The modified diffraction coefficients are 



n 2 cos(<5) 5 



£^(A£+e^)[A£+2e^ + Aecos(2g )] 2 
[Ae 2 +2Aee i +2£^+A £ (A£e i ) cos(20 o )] 2 

_ n 2 e sm($f _ 

e± [Ae+2e^+Ac cos{26 )} 

As 2 +2Aee±+2s' 2 L -r-Ae(Ae+2e±)cos{28( ) ) ' 



(5) 



D t =/= D x involves the absence of radially symmetric spa- 
tial solitons, with ellipticity (ratio between waists across 
t and x, respectively) given by Q = (Dt/Dx) 1 / 4 (see be- 
low). Noteworthy, when the birefringence Ae — > 0, it 
is D x = D t = 1 and isotropic propagation is retrieved. 
Figure [2] plots these quantities versus 8q for a highly- 
birefringent NLC.0] We need to stress that paraxiality 
in the walk-off system does not imply paraxiality in the 
original reference (x, y, z), as witnessed by the fact that, 
when re- writing equation i n (x, y, z) the second-order 
derivatives with respect to z re-appear. Since J3J holds 
for any walk-off, this treatment can model wide angle 
steering of spatial solitons, as e.g. obtainable by exploit- 
ing the voltage dependence of 5 in the NLC geometry. 

The molecular director n lies in the (y, z) plane and 
can be expressed in terms of angle = 8q + 'F (see fig- 
ure Using a multiple scale expansion, for the e-wave 
perturbation along t at the lowest order in 'F we have 
Ssu = i-8e-i= AeT(6 )^ with (see figure© 
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2e ± (Ae + e ± )sin(26» ) 



(Ae + e ± ) 2 +e 2 ± + [(Ae + £ ± ) 2 - e 2 ] cos(20 o ) 



(6) 



The NLC orientation is described by the functional F = 
Fk + Fjif + F opt , with Fk (Frank-term) accounting for 
the elastic properties of the NLC and Frf for its director 
distribution (and hence 9q) due to an externally applied 
(electric or magnetic) field in the absence of light. [l5Llr^ | 
For a dominant e-wave, the optical contribution to the 
energy F is F opt = ~[e ± \E e \ 2 + Ae(n • t) 2 |£ e | 2 ]/4. 

In the single-constant approximation (i.e., K = K± = 
K2 = K3 for molecular splay, bend and twist, respec- 
tively), from the Frechet derivative of F and at the 
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FIG. 2: (Top) Dimensionless anisotropy coefficients Vs 9o, 
see text. (Bottom) Walk-off angle 5 Vs #o; the insets show 
the soliton transverse intensity profiles for various 9o (labeled 



in degrees). Parameters: e±_ 
(max(n e ) — no + 0.5). 
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lowest-order in ^ we get: [l5llrRll^ | 

KVH - A(6o)V + ^ sin[2(0 o - 5)}\E e \ 2 = 0. (7) 

A(9q) is determined by cell geometry and bias. For the 
configuration in figure^ #o is determined by the voltage- 
driven elevation r] in the middle of the cell and by the 
azimuth £ due to NLC anchoring (at the interfaces defin- 
ing the cell) with respect to z [18|. We obtain 



A(e ) 



EqAsrf 

cos(£) 2 



sin(26» ) 
2#n 



cos(29 ) 



, (8) 



with Asrf the low-frequency (relative) permittivity, L 
an effective cell thickness over which the voltage V is 
applied and cos(# ) = cos(£) cos(r?o)- 
Spatial solitons. Optical spatial solitons or nematicons 
[l2j are solutions of equations (QJ) and in the form 
E e = (2Z /n e ) 1 / 2 U(x,t)exp(i(3s), with U 2 the intensity 
profile and (3 the "nonlinear wavevector," and d s ^> = 
(hereafter, we will always take 9 S 4' = in J7J), since 
the optical field is slowly varying along s, see also [lij). 
Self trapped beams travel along s, while their phase pro- 
file is orthogonal to the plane kon e z + (3s = [kon e + 
f3cos(5))z — (3sin(6)y ^constant, implying that for small 
(3 (low-power solitons) the phascfront corresponds to a 
plane wave propagating along z, gradually tilted towards 
s as the power increases. In other words, the nonlinearity 
tends to reshape the extraordinary wave into an ordinary- 
like configuration, by distorting the phase-fronts towards 
the plane orthogonal to the Poynting vector. The exact 
soliton profiles U can be obtained numerically. Neverthe- 
less, relevant insights can be obtained in the highly non- 
local limit, [H H3 as applicable to NLC. For NLC 
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FIG. 3: (Left) 5 Vs k and (Right) 8 Vs ellipticity, showing 
the soliton families attainable for the same walk-off. In the 
adopted units k measures the power in mW needed for a 1/im 
waist soliton. 



as in actual experiments, in fact, the perturbation \& ex- 
tends far from the excitation, so that the beam essentially 
experiences an index perturbation with a parabolic-like 
shape. Writing ip = ip + ^(x 2 +t 2 ), the equation for 
U can be analytically solved by separation of variables, 
yielding a wide class of self-trapped solutions including 
higher-order and breathing ones. The simplest profile is 
gaussian with intensity profile: 



P 2 r_P ( J_ 




(9) 



where P is the soliton power, and k is the constant of 



the existence curve: Pwq 



with wq the intensity 



(1/e) waist in the isotropic limit (D x = Dj = 1). It 
is k = 2Kn e c\ 2 ^/ D x D t /ir&E 2 Tsin(26 a - 25). In deriv- 
ing (O, we used V2 = — 2"oAe sin[2(0 — S)]/8Kn e c, as 
found from J2J) with lo the peak intensity when A — > 
(highly nonlocal regime). The self-trapped beam travels 
at any angle S(0q) with a gaussian profile and ellipticity 
Q = (Dt/Dj;) 1 / 4 , as anticipated. As in the case of MI (see 
below), two solitons (with different k) exist for the same 
5: their family is multivalued (one for each k, spanned 
by the power P) with respect to walk-off 8 (and unfolded 
by 6*o, i.e. propagating in different planes), as visible in 
figure |3 Noticeably, a straightforward generalization of 
the theory in [ilj (see also jj^) enables us to state that 
such solitons (i.e. those of eqs. (@J and Q) are uncondi- 
tionally stable. Indeed, the Hamiltonian for the system 
is written as H = J D t \d t E e \ 2 + D x \d x E e \ 2 dxdt + H nl 
where H n i the nonlinear nonlocal part which is identical 
to the isotropic case, and bounded from below as shown 
in 21] ; the remaining part is also bounded because, from 
Eq. QJ, D t and D x are not smaller than e±/(e± + Ae) 
(see e.g. figure |2J and this implies the boundness of H , 
and hence the stability. 

One- dimensional modulational instability. A sufficiently- 
wide elliptic beam propagating along z in the nonlinear 
sample approximates well a one-dimensional plane-wave. 
Owing to reduced diffraction across the major axis of 
the ellipse, in fact, the nonlinearity acts mainly in one- 
dimension. As confirmed by experiments, [l(l | MI causes 
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FIG. 4: (a) Maximum Mi-gain Vs 9 ; (b) as in (a) Vs walk- 
off (5; (c) maximally amplified spatial harmonic Vs 8o; (d) 
as in (c) Vs S. Thick (thin) lines refer to I = x (I = i). 
Parameters: V = IV, L = 75/im, £ - 
ex = 2.25, Ae = 1.75, E = 5 x 10 4 Vm" 
Eq.@. 



0, K = 10~ LL N, 
and A is given by 



the initially uniform beam profile to break-up of into a 
periodic pattern and, eventually, into periodically-spaced 
filaments. The instability can be intuitively described as 
the (selective in transverse spatial frequency) amplifica- 
tion of small amplitude noise superimposed to the input 
beam. For simplicity, we consider two limiting cases: an 
input ellipse with long axis oriented i) along x (I — x) 
or ii) along t (I — t). Plane- wave noise components of 
wavevector k x ,t will grow in amplitude along s with gain 
exp\gi(ki)s], being I — {x,t}: hence filaments form along 
the walk-off direction. A standard approach l23| pro- 
vides: 
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rDxh EikZe Ae 2 Tsm[2(9o-5)] 



2kon e cos(<5) 



2A 



2Kkf 



- Dikf 



(10) 

being Eq the peak amplitude of the pump (input) field. 
Ea. (|10|l is the generalized expression of Ml-gain in the 
presence of both non-locality and anisotropy, and reduces 
to the known result m the isotropic regime. The MP 
gain is spectrally affected by both the spatial orientation 
of the input ellipse (with major axis parallel to either x 
or t axes in cases i) and ii), respectively) and 9 . Fig- 
ure 0] shows the calculated gain profile for either orien- 
tations Vs #o (i-e. by varying the cell bias) and typical 
NLC parameters. Notably, the peak gain corresponds 
to the maximum walk-off; when graphed Vs 5, however, 
MI is multivalued: for the same S two distinct wave- 
patterns can emerge, belonging to different planes (t, z) 
and corresponding to different f3 . In an advanced stage 
of Mi-induced filamentation, this implies the possibility 
of angularly steering an entire array of regularly spaced 
(soliton) channel waveguides by acting on 9q. 

In conclusion, by developing a comprehensive model 
for nonlinear wave propagation in the presence of signifi- 
cant walk-off, we predict the existence of (multi-value) 



spatial solitons and modulational instability in highly 
birefringent and non-local media. Self-trapped beams 
travel at arbitrarily large walf-off angles, which deter- 
mine their elliptic intensity profile. The results hold valid 
for individual solitons and arrays of filaments as gener- 
ated through modulational instability, and can be read- 
ily extended to other media. In NPC, where walk-off 
can be adjusted by an external voltage, wide angle bias- 
controlled steering of ultra-thin anisotropic solitons could 
be effectively implemented, leading to applications such 
as optical information processing in space (e.g. spatial 
de-multiplexing) and optical tweezers. 
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